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Domelike magnetic-flux-density distributions previously have been observed experimentally and 
analyzed theoretically in superconducting films with edges, such as in strips and thin plates. Such 
flux domes have been explained as arising from a combination of strong geometric barriers and weak 
bulk pinning. In this paper we predict that, even in films with bulk pinning, flux domes also occur 
when vortices and antivortices are produced far from the film edges underneath current-carrying 
wires, coils, or permanent magnets placed above the film. Vortex-antivortex pairs penetrating 
through the film are generated when the magnetic field parallel to the surface exceeds H c \ + K c , 
where H c \ is the lower critical field and K c = j c d is the critical sheet-current density (the product of 
the bulk critical current density j c and the film thickness d). The vortices and antivortices move in 
opposite directions to locations where they join others to create separated vortex and antivortex flux 
domes. We consider a simple arrangement of a pair of current-carrying wires carrying current Iq in 
opposite directions and calculate the magnetic-field and current-density distributions as a function 
of Io both in the bulk-pinning-free case (K c = 0) and in the presence of bulk pinning, characterized 
by a field-independent critical sheet-current density (K c > 0). 

PACS numbers: 74.25.Sv,74.25.Nf,74.78.-w 



I. INTRODUCTION 

The hysteretic penetration of magnetic flux into a 
superconductor has long been a subject of experimen- 
tal and theoretical interest. One phenomenon that 
has often been observed in fiat samples in an increas- 
ing perpendicular applied magnetic field is a domelike 
magnetic-flux-density distribution centered in the mid- 
dle of the sample, surrounded by a superconducting flux- 
free zone. This phenomenon has been investigated ex- 
perimentally and interpreted theoretically in both type- 
1 superconductorsi*2iiMi£i£iL2iil and weak-pinning type-11 
superconductors, 10 : 11 : 12 ! 13 . 14 ! 15 . 16 ! 17 ! 18 ! 19 ! 20 In these pre- 
vious investigations, the flux domes have been explained 
as arising from an energy barrier of geometric origin at 
the edge of the sample. Once this barrier is overcome, 
flux tubes or vortices escape from the edge and are driven 
to the middle by Meissner screening currents, which flow 
on the sample's surface. 

In this paper we consider a type-II superconducting 
film without edges and predict that when it is subjected 
to local magnetic fields produced by current-carrying 
wires above the sample, vortex and antivortex magnetic- 
flux domes are produced in the film when the currents 
in the wires are sufficiently large. When bulk pinning is 
weak, we predict that the vortex and antivortex domes 
are separated from each other, but as the bulk pinning 
increases, the space between the vortex and antivortex 
domes shrinks to zero. 

To investigate these effects in an easily calculable ge- 
ometry, we consider in Sec.[n]a simple model in which the 
local magnetic fields are produced by a pair of infinitely 
long straight wires. The resulting two-dimensional ge- 



ometry allows us to calculate all the magnetic-field and 
sheet-current distributions analytically. In Sec. Ill Al we 
discuss the Meissner-state response of the film before any 
penetration of vortices into the film, and in Sec. IIIBl wc 
discuss the distributions produced by vortices and an- 
tivortices that have penetrated through the film thick- 
ness. We then present the magnetic-field and sheet- 
current distributions associated with the flux domes, 
both in the absence of bulk pinning [Sec. Ill C) and in 
the presence of weak [Sec. Ill D| and strong [Sec. Ill E| 
bulk pinning. In Sec. IIIII we summarize our results, dis- 
cuss the generality of the predicted effects in more easily 
realizable experimental configurations, consider similar 
phenomena in type-I superconductors, and discuss pos- 
sible extensions of this work. Calculations of screening 
effects in the Meissner state are presented in Appendix 
A, and derivations of the complex field and complex po- 
tential are given in Appendix B. 



II. SUPERCONDUCTING FILMS AND LINEAR 
WIRES 

We consider a simple geometry in which vortex and an- 
tivortex flux domes are produced far from the film edges. 
For simplicity we consider an infinite type-II supercon- 
ducting film and a pair of infinitely long current-carrying 
wires, as shown in Fig. 1. A wire, carrying current l\ = Iq 
parallel to the z axis, is situated at (x, y) = (0, y{), where 
yi > 0. The return current 1% = —Jo is carried by a sec- 
ond wire at (x, y) — (0, y%), where y% > y\. The wire 
radius r w is assumed to be much smaller than either 
y\ or the wire separation yi — y±. A superconducting 
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FIG. 1: Wire 1 at (x,y) = (0, j/i) carries current Ji = Jo in 
the z direction and wire 2 at (x, y) = (0, 3/2) carries current 
I2 = —Jo above an infinite type-II superconducting film in 
the xz plane. 



positive infinitessimal, since d <C 2/1,) the perpendicular 
and parallel magnetic fields H y (x,0) = ReH(x ±ie) and 
H x (x, ±e) = Im7Y(x ± ie) are obtained from Eq. as 



ff y (as,0) = #o a (M) + 



2tt 



fl»(a;,±e) = H 0x (x, 0) =F K z (x)/2, 



(4) 
(5) 



where Ho y (x,0) = Re H.q(x), Hq x (x,0) = Im7io(x), and 
P denotes the principal value integral. The complex po- 
tential is defined by Q(C) — /• W(C')^C'; an d the contour 
lines of the real part of correspond to the magnetic 
field lines of H{Q. 



A. Meissner-state response 



film, infinitely extended in the xz plane, is situated at 
—d/2 < y < d/2, where the film thickness d is larger 
than the London penetration depth A and d <C y\. As 
discussed in Appendix A, these conditions guarantee that 
when the film is in the Meissner state, the magnetic field 
below the film is negligibly small. 

Flux pinning in the film is characterized by the critical 
current density j c , which is assumed to be constant (in- 
dependent of magnetic field), as in Bean's critical state 
model, 2 - 1 and isotropic (independent of vortex direction) 
However, the relevant physical quantity here is the crit- 
ical sheet-current density K c — j c d. Since we are inter- 
ested in the case for which d <C yi, in the following we 
ignore the finite thickness of the film, as this simplifica- 
tion allows us to obtain simple analytic expressions for 
the magnetic-field and current-density distributions. 

We introduce the complex field TC(Q = H y (x,y) + 
iH x (x,y), which is an analytic function of £ = x + iy 
except for poles at £ = £1 = iyi and £ — C2 = iyi and a 
branch cut at y = 0. The Biot-Savart law for the complex 
field is given by 



00 c-« 



where 



Wo(C) = 



Jo 



1 



Jo 



1 



2tt C - iyi 2vr C - iyi ' 



(1) 



(2) 



is the complex field arising from the pair of wires alone 
(see Fig. 1) and K z (x) is the sheet current in the film. 
The complex potential describing the field generated by 
the wires alone, defined by 0o(C) = / C H (C)dC, is 



7T V I + iQ y 2 



(3) 



and the contour lines of the real part of Go(C) corre- 
spond to the magnetic field lines of Wq(0- At the up- 
per (C = x + ie) and lower (£ = x — ie) surfaces of the 
superconducting film, (where we take e = d/2 to be a 



We first consider the magnetic-field distribution when 
the film is in the (vortex-free) Meissner state and wires 
I and 2 carry dc currents I\ = Iq > and 1% = — Iq < 
after monotonically increasing in magnitude from zero. 
As discussed in Appendix A, when the current Iq is small 
and the thickness d is larger than A, the magnetic field 
is practically zero below the film, where y = fm£ < 0. 
The field distribution above the film can be obtained by 
adding to Hq(() and Go(0 the contributions Wj(C) and 
Qi(0 due to image wires at Q = — £i = —iyi and £ = 
— (2 = —iyi carrying currents — Iq and +Iq, respectively: 



Gi(0 = 



1 



r 



1 



2tt C + iyi 2tt( + iy 2 



In 



1 - i 



C/V2 \ 



(6) 
(7) 

for Im £ > 0, 

for Im ( < 0. 

(8) 

The subscript M is a reminder that this field describes 
the Meissner-state response to the applied field given in 
Eq. 0. Note that 



A-iC/yiJ 
The resulting complex field Hm(0 = ^o(C) 

.Io ( _ yi yi 

H M (0 







Jo 



H Mx {x,e) = — 



1)2 



7T \X A + 



vl 



1)1 



(9) 



As is usual using the method of images, Hmx{x,€) — 
2H Ox (x,0) = 2lmH (x) [see Eq. ©]. The corresponding 



complex potential Gm (C) 
given by 

.Jor 



jin M (C)dCior ImC>0is 



5m(0 



arctan 



( — ") — arctan (—) 
v yi ' V 2/2 / 



(10) 



We may take Gm{Q = for Im£ < 0. The perpendicular 
magnetic field and sheet-current density are thus given 
by Hm v {x,0) = and 



K Mz (x) = — -(- 

IT \C 



2/1 



1)2 



1)1 



(11) 
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FIG. 2: Contour plot of the real part of the complex poten- 
tial Qm{x + iy) vs x and y for y\ — 1 and j/2 = 2. The 
contours correspond to magnetic field lines of the complex 
field TLm(x 4- iy) describing the Meissner-state response of 
the superconducting film to currents in the wires shown in 
Fig. 1. The contours near the wires at £ = iyi and £ = ij/2 
are nearly circular. 




FIG. 3: Parallel magnetic field Hmx(x,£) / H c \ vs x (solid 
curve) in the Meissner state at the top surface of the super- 
conducting film when Io ~ I&, Hmx(0,€) = Hci, y\ = 1, 
and j/2 = 2. The dashed curve shows the corresponding 
sheet-current density Kmz{x)/H c i vs x. Note that Kmz{0) = 
—H c i. 



The net current induced in the superconducting film is 
I-x KMz(x)dx — 0, as expected. 

The maximum magnetic field parallel to the top surface 
of the film is 



H 



M. 



; (0,e) = ^(l-l), 

IT \Vl W>> 



(12) 



• 2/i 2/2 ' 

and the maximum magnitude of the sheet-current density 



\k m m\ = H--- 

7T Vyi y 2 



(13) 



We now define two critical currents as follows. Let 
I c i denote the value of Io at which Hm x (0, e) in Ea. ll2|) 
reaches the lower critical field H&, 



Id 



2/2-2/1 



(14) 



and let I c o denote the value of Io at which |-Kmz(0)| in 
Eq. i|13|) reaches the critical sheet-current density K c , 



IcO 



7T^c2/l2/2 

2/2-2/1 



(15) 



Both of these critical currents play important roles in 
determining the details of how magnetic flux penetrates 
through the film. In the following sections we discuss in 
turn the cases for which K c = (I c o = 0) and K c > 
{Ico > 0). 



B. Vortex-generated fields and currents 

Consider a closely spaced row of vortices along the z 
axis carrying magnetic flux $' per unit length up through 
the film. Ignoring spatial variation on the scale of the 
London penetration depth A or the intervortex spac- 
ing, the magnetic field in the space y > (y < 0) ap- 
pears as if produced by a line of positive (negative) mag- 
netic monopoles. At a distance r from the origin, the 
magnitude of the magnetic field is h = <£>'//io7rr. Ex- 
pressing this result in terms of a complex magnetic field 
^u(C) — H V y(x,y) + iH vx (x,y) and extending it to a 
distribution of vortices or antivortices generating a mag- 
netic field H vy (x,0) in the plane of the film, we see that 
the vortex-generated complex magnetic field can be ex- 
pressed as 



n v (0 = ±~ J 



+°° H vy (x,0) 
du—? . 



(16) 



where £ = x + iy and the upper (lower) sign in Eq. i|lu|) 
holds for y > (y < 0). The corresponding vortex- 
generated sheet-current density K vz (x) — H vx (x — ie) — 
H vx (x + ie) is 



K vz (x) 



du 



2H vy (x,0) 



(17) 



while the Biot-Savart law yields another relation between 
K vz (x) and H vy (x,0): 



2H vy (x,0) = -[ 



+°° , K vz (u) 

du . 

x — u 



(18) 



In the following sections, the complex field can al- 
ways be regarded as a linear superposition of the 
Meissner-state and vortex-generated complex fields: 
U(0 = H M {Q +%u(C). Th e complex potential Q(Q = 
jf e H(C)dC can be written as G(() = G M (C) + G v (0- 
Similarly, the sheet-current density can always be ex- 
pressed as K z (x) = Kmz{x) + K vz (x). 



C. Flux domes in the absence of bulk pinning 

In bulk-pinning- free films (K c — 0), the first vortex en- 
ters the film at x = 0, where the maximum magnetic field 
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at the top surface is equal to the lower critical field H c i 22 
This occurs at the current Iq = J c i, given in Eq. i|14|) 
(see Fig. 3). An initially tiny vortex loop expands in ra- 
dius, and a portion of the loop is driven to the bottom 
of the film surface, where it annihilates, resulting in a 
separated vortex-antivortex pair. The vortex (antivor- 
tex) carries magnetic flux </>o in the +y (— y) direction, 
where <pQ = h/2e is the superconducting flux quantum. 
Responding to the Lorentz force Kmz( x )<Po, the vortex 
moves in the x direction until it comes to rest at the 
point x — x — where Km z {x) = 0, as can be 

seen from Eq. and Fig. 3. The antivortex moves 
in the opposite direction and comes to rest at the point 
x = —xq. 

For increasing values of Iq in the range I c \ < Iq < 
2I C \ , the magnetic field distribution perpendicular to the 
film can be characterized as having a positive vortex- 
generated magnetic flux dome in the region b < x < a, 
where < b < xq < a, and a negative antivortex- 
generated flux dome in the region —a < x < —b. The 
complex magnetic field H(() — Hm{() +Tt v {() is given 
in Eq. (|B13|I but with K c = 0. Subtracting the Meissner- 
state complex field H.m((), we obtain the following ex- 
pression for the vortex-generated complex magnetic field, 




FIG. 4: Plot of the vortex dome's left boundary b vs Io/I c i for 
the bulk-pinning-free case, where b is in units of xo — <JyiV2, 
for Jci < Io < 2ici and y x /y 2 = 0.01, 0.03, 0.1, 0.3, and 1 
(bottom to top). The right boundary of the vortex dome is 
at a = Xo/b. 

{Io _ y 2 [l-cj>(x)/s 2 } \ 

2tt\ x 2 + y\ x 2 + y 2 )' 

\x\ < b or \x\ > a, " ( 23 ) 
-^(i^-5^).6<N<«. 



Z7T t 



Io (yi[Ti- 



<f>(0/ si] 2/ 2 [Ti + <KC)/s 2 ] 



C 2 + y 2 2 



}, (19) 



where Sj — ■J (a 2 + y 2 )(b 2 + y 2 ), and the upper (lower) 

sign holds when £ = x + iy is in the upper (lower) 
half plane. Writing </>(£) = <f>i(x,y) + i4> 2 {x,y), we find 
that 4>\(x,y) is an odd function of x and an even func- 
tion of y [(j>i(x,y) = -cj>x(-x,y) = <j>i(x,-y)], while 
(j) 2 (x,y) is an even function of x and an odd function of 
y [4>2(x,y) = (j) 2 (-x,y) = -<j> 2 {x,-y)). Thus H vy (x,y) 
is an odd function of x but an even function of y, while 
Hvx{x, y) is an odd function of y but an even function of 
x. Just above (below) the real axis, 



I (j>(x), b < \x\ < a, 



where 



V(a 2 ~x 2 )(b 2 -x 2 ), 



\x\ < b, 

4>(x) — ^ sgn(x)^(a 2 — x 2 )(x 2 — 6 2 ), b < \x\ < a, 
- \J (x 2 — a 2 )(x 2 ~ 6 2 ), \x\ > a. 

(21) 

From Eq. ifH^I we obtain the following values of H vy (x, 0) , 

H vx (x,e) = -H vx (x,-e), and K vz (x) = H vx (x,-e) - 



H. 



0, \x\ < b or |a;| > a, 

(x,0) = <| ^(^(JVa?) ( 22 ) 



io fyi[i-4>( x )/ s i] y2[i-<j>{x)/s 2 ] 



K vz (x) 



-( 

7T V 



x 2 + y\ 



x 2 + y 2 



)■ 



in( Vl 

. 7r \x 2 +yj 



< b or \x\ > a, 



(24) 



b < \x\ < a, 



As discussed in Appendix B, the requirement that 
j+00 j^ z ^ x ^ x _ q } eac i s to the condition that ab = Xq = 
2/12/2 when there is no bulk pinning. A second condition 
relating a , b and Iq follows from the requirement that 
H x (0, e) = H Mx (0, e) + H vx (0, e) = H cU which yields 



Io 



1 + 2/iy 2 (si + s 2 )/sis 2 



Id- 



(25) 



Combining these two conditions and eliminating a, we 
obtain the following connection between Iq and b, 



l + b( yi + y 2 )/V(yl + b 2 )(y 2 2 +b 2 ) 



)ld- (26) 



When Iq = I c i, a = b = xq = ^Jy\y 2 . As Iq increases 
above I c \ , b decreases monotonically, as shown in Fig. 01 
and the lower solid curve in Fig. ^2 I n the limit as 
Iq — > 2/ c i, we see that b — > 0, such that a = y±y 2 /b — > 00. 
For Jo > 2/ c i, the magnetic-flux-filled region extends 
from — 00 to +00, K z (x) = 0, and H(C) = 7~Co(() every- 
where, since the superconducting film is then completely 
incapable of screening the magnetic field produced by the 
two wires. 

Figures 5 and 6 show plots of H y (x, 0) = H vy (x, 0) and 
K z (x) vs x for several values of Iq in the range I c \ < Iq < 
2/ c i. 
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FIG. 5: The perpendicular magnetic 
H vy (x,0)/ H c i vs a; exhibits vortex and 
at b < | a; | < a in the absence of bulk 
1, 7/2 = 2, and x = ^/j/Iyi = 1.414, 
tiny domes are centered at ±xo- As Io 
become taller and wider, as shown for 
6=1 (dotted), Jo = 1.2127 c i when b = 
Jo = 27 c i when 6 = and a = oo (long 
the superconducting film is no longer 
K z (x) — everywhere, and H y (x,0) — 



field H y (x,0)/H cl = 
antivortex flux domes 
pinning. With y\ = 
when J = 1.0017 c i, 
increases, the domes 
Jo = 1.026 Jci when 
0.5 (short dash), and 
dash). For J > 2J c i, 
capable of screening, 
#02,(2;, 0) [Eq. ©]. 
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FIG. 6: Sheet-current-density K z (x) / H c \ vs x in the absence 
of bulk pinning for yi = 1, j/ 2 = 2, and so = y/yiy-i = 1.414. 
When Jo = J c i, 7T z (a;) (solid) is equal to Kmz(x), shown in 
Fig. 3. With increasing Jo, regions of K z (x) = develop 
underneath the vortex and antivortex domes at b < \x\ < a, 
beginning at x = ±xo- Shown are results for Jo = 1.026J c i 
when 6=1 (dotted) and Jo = 1.212J c i when b — 0.5 (dashed). 
For Jo > 2 Jci, the superconducting film is no longer capable 
of screening, and K z (x) — everywhere. 



Figure 7 shows plots of H x (0,±e) = Hmx{0, ±e) + 
H vx (0,±e) vs Jo- For < Jo < 7 c i, we have 
H x (0,e) = H Mx (0,e) - H cl I /I cl and H x (0,-e) = 
Hjrfx{0, — e) = 0. For I c \ < I < 27 c i, we have 
H x {0, e) = H Mx {0, e) + H vx (0, e) = H cl and H x (0, -e) = 
H vx (0,-e) = H cl (I /I cl - 1), where ff Mx (0,e) = 
Hcilo/ Id and Jr ox (0,±e) = ±H cl (l -J /I d ). When 
I = 27 cl , H x {0,±e) = H cl . For J > 7 cl , 77 x (0,±e) = 
7?oa:(0,0) = 77 c i7o/27 c i, and in addition H x (x, ±e) = 
HMx{x,±e) + H vx (x,±e) — Hq x (x,0) for all X. The 
macroscopic magnetic-field distribution is then the same 
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FIG. 7: Parallel magnetic fields at a; = at the top and 
bottom surfaces H x (0,e) and H x (0, — e) vs Jo in the absence 
of bulk pinning. When < 7 < 7 c i, H x (0,e) = H M x(0,e) 
[Eq. ©] and 77 x (0, -e) = H Mx {Q, -e) = 0. When 7 cl < 
Jo < 27d, JJx(0,e) = JJ A /x(0,e) +77 ra (0,e) = H cl , and 
H x (0,-e) = H vx (0,-e) = -H vx (0, e). When J > 2J cl , 
ff„(0,±e) = JJox(0,0) = H cl I /2I cl [Eq. ©]. 



-K (0) 




FIG. 8: The magnitude of the sheet-current-density at the 
origin —K z (0) = H x (0, e) — H x (0, — e) vs Jo in the absence of 
bulk pinning (see Fig. 7). 



as it would be if the superconducting film were absent. 

Figure 8 shows a plot of —K z (0) vs Jo, where K z (0) = 
H x (0,-e) - H x (0,e). For < 7 < 7 cl , -K z (0) = 
Hdlo/Ici, and for J cl < I < 27 cl , -K g (0) = 77 cl (2 - 
lo/Ici)- For Jo > 2/ c i, K z (x) — for all x, i.e., every- 
where in the film. 

The vortex-generated complex potential Gv(C) = 
G{() - Sm(C) [see Eqs. JTHJl and QETty ] in the absence 
of bulk pinning is 



Qv(0 



h 

2tt 



\9v(£,Vi) - 9v(t,Va)] 



(27) 



where 

9v{(,y) 



=pi arctan(C/y) 



T [a 2 y 2 E(9 7 k) + y 2 (b 2 + y 2 )F(9, k) 

asy 

-(a 2 + y 2 )(b 2 + y 2 )n(9, -b 2 /y 2 ,k)}, (28) 



arcsin(£/&), and 
b/a, 



(29) 
(30) 
(31) 



where E, K, and II are incomplete elliptic integrals and 
the upper (lower) signs hold in the upper (lower) half £ 
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FIG. 9: Contour plot of the real part of the vortex- generated 
complex potential in the absence of bulk pinning [Eq. 1271 1 
Q v (x + iy) vs x and y for y\ = 1, j/2 = 2, a = 2, and 6=1, 
as for the dotted curves in Figs. 5 and 6. The contours corre- 
spond to magnetic field lines of the vortex-generated complex 
field H v (x + iy) given in Eq. IT9l . 



plane. Shown in Fig. 9 is a contour plot of the real part 
of G v {x + iy). These contours correspond to the mag- 
netic field lines of the vortex-generated magnetic field. 
The magnetic field flows in a generally counterclockwise 
direction, carried by vortices in the region b < x < a 
up through the film and by antivortices in the region 
— a < x < — b back down through the film. A contour 
plot of the real part of G(x + iy) = Gm (x + iy)+Gv(x + iy) 
would also show the magnetic field lines generated by the 
two wires, as in Fig. 2. 



D. Flux domes in the presence of weak bulk 
pinning 

In superconducting films in which bulk pinning is 
present and is characterized by a field-independent crit- 
ical sheet-current density K c — j c d > 0, the first vortex 
enters the film at x = 0, when the maximum parallel 
magnetic field at the top surface fjjv/x(0, e) [Eq. iQ ] is 
equal to the lower critical field H c iml As in Sec. Ill Cl this 
again occurs at the current Iq — I c \, given in Eq. (|14f> . 
According to critical-state theory^24 this vortex ad- 
vances toward the bottom surface at the leading edge 
of a curving flux front of nearly parallel vortices of thick- 
ness dp = [Hmx(0, c) — Hci]/jc, and it reaches the bottom 
surface when d p — d, i.e., when Hmx(0,c) = H c \ + K c 
or I = I cl + J c0 [see Eqs. lfT2 ]) . and JTSJl]. The 

positive end of the first vortex is then at x = x p and its 
negative end is at x = —x p , where x p is the solution of 



HMx(x p ,e) = H c i [see Eq. JjjJ], It can be shown that 
x p ^> d when d <C y\ except when K c <C H c \ . 

Once the first vortex reaches the bottom surface, the 
portion at x ~ annihilates with its image, and the 
vortex divides into two halves. The half in the region 
x > 0, which we call a vortex, carries magnetic flux (f>Q 
up from the bottom to the top surface, and the half in 
the region x < 0, which we call an antivortex, carries 
magnetic flux 4>$ down from the top to the bottom sur- 
face. Since Hmx(0, — e) = and Hm x (0, e) = Hci + K c 
at Iq = I c i + I c o, the sheet-current density is initially 
Km z (0) = —Hci — K c - Because |ivMz(0)| > K c , the vor- 
tex separates from the rest of the flux front and is driven 
in the x direction by the Lorentz force Kmz{x)4>o until it 
comes to rest at x = x c , where \Kmz{x c )\ = K c [see Eq. 
Ullfl ] and the Lorentz force is balanced by the pinning 
force. Similarly, the antivortex moves in the opposite di- 
rection and comes to rest at the point x = —x c . As we 
will show below, < x c < xq = y/yiy2 when K c > [see 
Fig. 10]. 

For increasing values of Iq in the range I c i + I c o < h < 
2/ c i+/ c o, the magnetic field distribution perpendicular to 
the film can be characterized as having a positive vortex- 
generated magnetic flux dome in the region b < x < 
a, where < b < x c < a, and a negative antivortex- 
generated flux dome in the region —a < x < —b. The 
complex magnetic field H(() — 7Ym(C) + ^«(C) is given 
by Eq. (IJj 13|) . Subtracting the Meissner-state complex 
field 7Ym(C)j we obtain the following expression for the 
vortex-generated complex magnetic field, 



«-(0 = £{ 



»l[T* + ^(C)/»l] y 2 [T« + 0(C)/S2] 



C 2 + yl 



C + yi 



±iK c /2, 



(32) 



where Sj = ^JJa? 



¥j){b 2 + Vj), and the upper (lower) 

sign holds when £ = x + iy is in the upper (lower) half 
plane. Following a procedure similar to that used in 
Sec. Ill CI we obtain the following values of H vy (x,0), 
H vx (x,e) 
H vx (x, e), 



-H vx (x,-e), and K vz (x) = H vx (x,-e) 



H vy (x,0) 



In ( Vi 

2tt \^si(:e 2 - 



< b or Id > 



+vD 



(33) 



V2 



S2(x 2 +vi) 



4>(x), b < \x\ < 



Iq fyi[l-4>(x)/si] y 2 [l - (j)(x) / s 2 ] 



H vx (x,e) = < 



-Kc/2, 



„2 + y 2 



X- 



— ( m 

7T I x 2 +y'f 



\x\ < b or |x| > a, 
-V- 1 + K c /2, 



(34) 



b < \x\ < a, 
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FIG. 10: Plot of x c = Xc/xo vs K c jH c \ in the presence of bulk 
pinning, where x c is the position of the first entering vortex 
when Jo just exceeds J c i + J c o and xo = ^/j/ij/2- Curves are 
shown for j/i/j/2 = 0.01, 0.03, 0.1, 0.3, and 1 (bottom to top). 



h /yi[l-<KaO/si] y 2 [l-(j)(x)/s 2 



K vz {x) 



-( 

7T V 



^2 



^2 _|_ ^2 

x| < b or |x| > a, 

(35) 

- K c , b< \x\ < a, 



X -{-yf x 2j ty 

For I c \ + 7 c o < Io < 2I c i + J c o , the requirement that 
f*™ K z (x)dx — leads to Eq. 1B24(> . A second condi- 
tion on a, 6, and Jo follows from the requirement that 
H x (0, e) = H cl + K c , which from Eq. yields: 



H x (0,e) 



Io_ \ (1 + ab/si) (l + ab/s 2 ) 
2tt 



Hci 



Vi 

K r . 



!J2 



Kc 
2 



(36) 



Elimination of Jo between Eqs. (IB24II and l|36(l yields the 
equation 



= K c /2H cl 

{V2 - yi)sis 2 + ab(y 2 s 2 - yisi) 1 + K c j1H c \ 



y\y-i{y\s 2 - y 2 si) 



(37) 



Numerical solutions of Eqs. I|3()|) and l|37|) yield a and b 
as a function of Jo for any given value of K c . As dis- 
cussed above, when J just exceeds J c i + J c o, the first 
vortex (antivortex) comes to rest at x c (-x c )- The equa- 
tion determining the value of x c = x c /xo — x c / ' ^Jy\y 2 
can be obtained from Eq. i|B24(l by setting a = b = x c 
and making use of Eqs. (|14fl and l|15|) : 



(1 



K c /H c 



(x 2 c +yi/y 2 )(x 2 c +y 2 / yi ) 1 + K C /H C 



(38) 



Figure 10 shows plots of x c vs K c / H c \ for various values 
of yi/y 2 , obtained by numerically solving Eq. ll3*5|) . Note 
that for each case x c = 1 when K c = 0, and x c — > when 
K c /H cl -» oo. For Jy c /JJ c i > 1, 



So « i/V(i + yi/y2 + w/vi)(Kc/H cl ). (39) 

Numerical solutions for a and b when J cl + J c0 < Jo < 
2J c i + J c o are shown by the dotted curves in Fig. 11 for 



FIG. 11: Plots of a (upper curves), the right boundary of the 
vortex dome, and b (lower curves), the left boundary, vs Jo/Jci 
for yi — 1 and y2 = 2: K c — (solid curves separated by the 
dot at a = b — x c — 1.414), K c = JJ c i/2 (dotted curves and 
dot at x c — 0.827), K c — JJ c i (short-dashed curves and dot 
at x c = 0.651), and K c = 2JJ c i (long-dashed curves and dot 
at x c = 0.493). Note that b = when J /J c i > 2 + K c /H c i. 
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FIG. 12: When Jo > J c i + J c o, the perpendicular magnetic 
field H y (x, 0)/ Hci = H vy (x, 0)/H c i vs x initially exhibits sep- 
arated vortex and antivortex flux domes at b < \x\ < a even in 
the presence of bulk pinning, as shown here for K c = H c i/2, 
yi = 1, ?/2 = 2, and x c = 0.827. When Jo is just above 
Jci + Jco = 1.5J c i, tiny domes are centered at ±x c . As 
Jo increases, the domes become taller and wider, as shown 
for J = 1.545J c i when b = 0.600 and a = 1.077 (dotted 
curves), Jo = 1.683J c i when b = 0.400 and a = 1.337 (short- 
dashed curves), Jo = 1.966J c i when b — 0.200 and a — 1.630 
(medium-dashed curves), and Jo = 2J c i + J c o = 2.5J c i when 
6 = and a = 1.933 (long-dashed curves). 



the case that y\ = 1, y 2 = 2, xo — V2, and K c = H c ±/2, 
such that J c o = J c i/2. Note that a = b = x c = 0.585xo = 
0.827 at Jo = J c i + J c o when vortex penetration first 
occurs. From Eqs. ffity . 1(14^ . and (fT5|l . we see that b = 
when Jo = 2J c i + J c o. When Jo > 2J c i + J c o, b remains 
equal to zero, and the value of a must be obtained from 
Eq. I|B25|) . Figure 11 also shows similar plots of a and b 
for K c = H c i and K c — 2H C \ 

Figures 12 and 13 show plots of H y (x,0) — H vy (x,0) 
and K z (x) vs x for several values of Jo in the range J c i + 
Icq < Io < 2J c i + J c q. 
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FIG. 13: Sheet-current density K z (x,Q) / H c \ vs x in the pres- 
ence of bulk pinning, shown here for K c = H c i/2, y± — 1, 
2/2 = 2, and x c = 0.827. When 7 = I c i + I c o, K z (x) (solid) 
is equal to Kmz(x), shown in Fig. 3. With increasing To, re- 
gions of K z (x) — —K c develop underneath the vortex and 
antivortex domes at b < \x\ < a, beginning at x = ±x c , 
while \K z (x)\ > K c in the vortex-free region —6 < x < b. 
Shown are results for To = 1.545/d when b — 0.600 and 
a = 1.077 (dotted curves), T = 1.683T c i when b = 0.400 
and a — 1.337 (short-dashed curves), To = 1.966T c i when 
b = 0.200 and a — 1.630 (medium-dashed curves), and 
T = 27 c i + T c0 = 2.500T c i when b = and a = 1.933 (long- 
dashed curves). For To > 2T c i + T c o, K z (x) = — K c in the 
region —a < x < a and |Tf z (:r)| < K c outside this region. 





H x (0,e) / 










/tip,-*) 



FIG. 14: Parallel magnetic fields at x = at the top and 
bottom surfaces H x (0,e) and H x (0,—e) vs To in the pres- 
ence of bulk pinning, shown here for K c = TT c i/2, such 
that T c i + T c o = 1.5T c i and 2T c i + T c o = 2.5T c i. When 
< To < Tei + Too, 7^(0, e) = H M x(0,e) [Eq. Q] and 
77^(0, -e) = H M *(0,-e) = 0. When I cl + T c0 < T < 
2T c i + T c0 , H x (0, e) = H Mx (0, e) + H vx (0, e) = H cl + K c , and 
H x (0, -e) = H vx (0, -e) = -H vx (0, e). When T > 2T ci + T cU , 
H x (0,±e)=H 0x (0,0)±K c /2. 



Figure 14 shows plots of H x (0,±e) = H Mx (0,±e) + 
H vx (0,±e) vs Jo- For < Iq < I c i + I c o, we have 
H x (0,e) = H Mx (0,e) = H cl I /I cl and 77,(0, -e) = 
H Mx (0,-e) = 0. For I cl + I c0 < I < 2J cl + 7 c0 , we 
have 77,(0,6) = H Mx (Q,e) + H vx (0,e) = H cl + K c and 
77,(0,-6) = H vx (0, -e), where 77 Mx (0,e) = H cl I /I c i 
and H vx (0,±e) = tH c1 (I q - I cl - I c0 )/I cl . When 7 = 
27 cl + 7 c0 , 77,(0, e) - 77 cl + 7^ c and 77,(0, -e) = H cl , 
For 7 > 27 cl + 7 c0 , 77,(0, ±e) = 77 ,(0,0) ± K c /2 = 




FIG. 15: Magnitude of the sheet-current density at the origin 
—K z (0) = H x (0, e) — H x (0, — e) vs To in the presence of weak 
bulk pinning, shown here for K c = H c i/2, such that T C 2 = 
2T c i - T c0 = 1.5T c i (see Fig. 14). 

77 cl 7 /27 cl ±7f c /2. 

Figure 15 shows a plot of —7^(0) vs 7o, where 7^2(0) = 
77,(0, -e) - 77,(0, e). For < 7 < 7 cl + 7 c0 , -K Z (Q) = 
(H cl +K c )I /(I cl +I c0 ), andfor7 c i+7 c0 < 7 < 2I cl +I c0 , 
-K z (0) = K c + H cl (2I cl +I c0 -I )/Ici- For 7 > 27 cl + 
ico, ~K Z {0) = K c . 

The vortex-generated complex potential ^„(C) = 
Q(C) ~ Gm(0 [see Eqs. JTOJ and (|B14|) ] in the pres- 
ence of weak bulk pinning (0 < K c < 77 c i) is, when 
I c i < 7 < 7 c2 , such that < b < a, 

Qv(C) = ^ [9v(C vi) - 9v(C, V2)] ± C, (40) 

where g v (CiV) is given in Eq. I)28|). Shown in Fig. 16 
is a contour plot of the real part of Q v (x + iy). These 
contours correspond to the magnetic field lines of the 
vortex-generated magnetic field. The magnetic field flows 
in a generally counterclockwise direction, carried by vor- 
tices in the region b < x < a up through the film and 
by antivortices in the region —a < x < —b back down 
through the film. A contour plot of the real part of 
Q(x + iy) = Gm{x + iy) + Q v (x + iy) would also show 
the magnetic field lines generated by the two wires, as in 
Fig. 2. 

When 7o > 27 c i + 7 c o, the gap between the vortex 
dome and the antivortex dome is closed (b — 0), and 
the magnetic-field distribution thus can be characterized 
as a dome of vortices in the region < x < a carrying 
magnetic flux up through the film and an adjacent dome 
of antivortices in the region — a < x < carrying an equal 
amount of magnetic flux back down through the film. 
The outer boundaries of these domes (±a) depend upon 
the values of Iq and K Cl and the magnetic-field and sheet- 
current-density distributions can be calculated as follows. 
The complex magnetic field H(() = Hm(() + 7~tv(() is 
as given in Eq. (|B13I) . except that 6 = and (/)((), S\, 
and S2 are now given by Eqs. (|B19l) . fB20j, and |B2lJ. 
Similarly, the vortex-generated complex magnetic field 
H V (C) is as given in Eq. and H vy (x, 0), H vx (x, e) = 
-H vx (x,-e), and K vz (x) = H vx (x 1 -e) - H vx (x,e) are 
as given in Eqs. (|34|) . and l|35|l . except that now just 
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FIG. 16: Contour plot of the real part of the vortex-generated 
complex potential in the presence of weak bulk pinning 
[Eq. I|4()p] Q v (x + iy) vs x and y for y\ = 1, 2/2 = 2, 
A" c = H c i/2, Jo = 1.683/ci, a = 1.337, and b = 0.4, as for the 
short-dashed curves in Figs. 12 and 13. The contours corre- 
spond to magnetic field lines of the vortex-generated complex 
field TC v (x + iy) given in Eq. 1321 . 



above or below the real axis, 



4>(x ± ie) 




< |x| < a, 
\x\ > a, 



(41) 




FIG. 17: When I > 2I cl + I c0 in the presence of bulk 
pinning, the perpendicular magnetic field H y (x,0)/H c i = 
H vy (x,0) / H c i vs x exhibits adjacent vortex and antivortex 
flux domes at < |x| < a, as shown here for K c = H c i/2, 
Ico = 0.57 c i, yi = 1, and 1/2 = 2. As Iq increases, the 
domes become taller and wider, as shown for Io = 2.5/ c i when 
a = 1.933 (solid curves), Io = 4.25J c i when a — 2.504 (short- 
dashed curves), and Io = 67 c i when a — 2.913 (long-dashed 
curves) . 




where 



fa) 



'y a 2 — 
-\x\ vr 



< \x\ < a, 
\x\ > a. 



(42) 



For given values of Iq and K c the value of a in the above 
equations is determined by Eq. (|B25|) , which follows from 



the requirement that f 00 K z (x)dx = 0, and by Eq. 136(1 
(but with 6 = 0), which follows from the requirement that 
H x (0, e) = H c i + K c . Numerical solutions for a obtained 
in this way are shown in Fig. llll vs Iq/I c i for I > 2/ cl + 
Ico (where 6 = 0) for the case of y± = 1 and yi = 2. 
Figures UTI and IT51 show plots of H y (x,0) — H vy (x,0) 
and K z (x) vs x for several values of Io > 2/ c i + I c q when 
K c = H cl /2. 

When 7 > 2I c i + J c o, such that a > and 6 = 0, 
the vortex-generated complex potential Gv(C) — G(0 ~ 
Gm{Q [see Eqs. iJTDJ, (|BF4) . and l|B22)l ] in the presence 
of bulk pinning is 



G v (0 



Jo 
2tt 



[g v (<;,yi)-g v ((,y2)}±^(, (43) 



FIG. 18: Sheet-current density K z (x, 0)/H c i vs x when Io > 
2/ c i + / c o in the presence of bulk pinning, shown here for K c = 
H c i/2, Ico — 0.5/ c i, yi = 1 and t/2 = 2. Shown are results for 
Io — 2.5/ci when a — 1.933 (solid curve), Jo = 4.25/ci when 
a = 2.504 (short-dashed curve), and Io = 6/ c i when a = 2.913 
(long-dashed curve), corresponding to the cases shown in Fig. 

El 



where 



9v((,y) = =R arctan(C/y) 
— arctanh 



a 2 -C 2 _ 
a 2 +y 2 

arctanh 



^/ a 2 +y 2 



(44) 



and the upper (lower) signs hold in the upper (lower) half 
£ plane. Shown in Fig. El is a contour plot of the real 
part of Gv(x+iy). These contours correspond to the mag- 
netic field lines of the vortex-generated magnetic field. 
The magnetic field flows in a generally counterclockwise 
direction, carried by vortices in the region < x < a 
up through the film and by antivortices in the region 
—a < x < back down through the film. A contour plot 
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FIG. 19: Contour plot of the real part of the vortex-generated 
complex potential in the presence of bulk pinning [Eq. (1431 1 
Q v (x + iy) vs x and y for yi = 1, y 2 = 2, K c = H c i/2, 
Jo — 2/ c i + J c o = 2.5J c i, a — 1.933, and b — 0, as for the 
solid curves in Figs. 1171 and 1181 The contours correspond 
to magnetic field lines of the vortex-generated complex field 
Hv{x + iy) given by Eq. 1521 . 

of the real part of Q(x + iy) = Gm{x + iy) + Q v (x + iy) 
would also show the magnetic field lines generated by the 
two wires, as in Fig. 2. 



E. Flux domes in the presence of strong bulk 
pinning 

We consider here briefly the case of superconducting 
films in which bulk pinning, characterized by a field- 
independent critical sheet-current density K c = j c d > 0, 
is so strong that K c ^> H c \. In principle, the process 
of vortex entry is qualitatively the same as discussed 
in Sec. Ill Dl However, in the limit that H c i/K c — > 0, 
the width of the current region J c i + I c q to 2J c i + I c q 
(in which there is a gap of width 26 between the vortex 
and antivortex domes) shrinks to zero, and I c q becomes 
the only critical current of practical interest. Essentially, 
as soon as Jo exceeds J c o, the vortices penetrating from 
the top surface divide in such a way as to produce ad- 
jacent vortex and antivortex flux domes in the regions 
< x < a and — a < x < 0, and the sheet-current den- 
sity is K z (x) = —K c in these regions. The perpendicular 
magnetic field H y (x, 0), sheet-current density K z (x) and 
vortex-generated complex potential Gv(x + iy) can be cal- 
culated as discussed in Sec. Ill Dl for the case that 6 = 0, 
and plots of all these quantities look very similar to those 
in Figs. E| HHI andCa 

Figure 20 shows plots of H x (0,±e) = H Mx (0,±e) + 



FIG. 20: Parallel magnetic fields at x = at the top and 
bottom surfaces H X (Q, e) and H x (0, — e) vs Jo for strong bulk 
pinning in the limit as H cl /K c — > 0. When < Iq < I c o, 
fl»(0, e) = H Mx (0, e) [Eq. jgj] and H x (0, -e) = H Mx (0, -e) = 
0. When 7 > I c0 , H x (0,±e) = H Ox (0, 0)± K c /2 [see Eq. ©]. 



-K z (0) 




FIG. 21: Magnitude of the sheet-current density at the origin 
— K z (0) = H x (0, e) — H X (Q, — e) vs Jo for strong bulk pinning 
in the limit as JJ c i /K c — ► (see Fig. 20) . 



H vx (0,±e) vs lo in the limit H cl /K c — * 0. For 
< Jo < Jco, we have H x (0,e) = H Mx (0,e) = K c I /I c0 
and H x (0,-e) = H Mx {0,-e) = 0. For J > J c0 , 
H x {0, ±e) = JJ 0k (0, 0) ± K c /2 = K c I /2I c o ± K c /2. Fig- 
ure 21 shows a plot of — K z (0) vs Jo in the same limit, 
where K z (0) = JJ x (0, -e) - H x (0,e). For < J < J c0 , 
-K z (0) = KJo/Iao, and for J > J c0 , -K z {0) = K c . 



III. SUMMARY AND DISCUSSION 

In this paper we have predicted that separated vor- 
tex and antivortex flux domes can be produced in weak- 
pinning type-II superconducting films subjected to lo- 
cal magnetic fields generated by current-carrying wires 
above the film's surface and far from the edges. To cal- 
culate these effects analytically, we have chosen an ideal- 
ized geometry of two parallel infinitely long wires above 
an infinite superconducting film. However, the basic phe- 
nomenon of the creation of separated vortex and antivor- 
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tex flux domes without nucleation at the film's edges is 
far more general than in the geometry we have considered 
here. 

For example, consider the case of a bulk-pinning-free 
superconducting film of finite size in the xz plane and a 
small coil at (x,y, z) = (0,yoiO) a short distance above 
the film. Suppose the coil produces a magnetic dipolc 
moment m = —im. When m is small, the film remains 
in the Meissner state, and if the magnetic field below 
the film is very small, the coil-generated dipole magnetic 
field and its image produce a magnetic field at the film's 
top surface Hm = xHm x + zHmz, where Hmx(x, e, z) — 
2m(j/o - 2a; 2 + z 2 )/Rq, H Mz (x,e,z) = ~6itixz/Rq, and 
Ro = \Jx 2 + + z 2 . The induced sheet-current den- 
sity is K M = xK Mx + zK Mz , where K Mx {x,z) = 
Hmz(x,€, z) and K~Mz(x,z) — —Hmx(x,€,z), such that 
V • Km = 0. In the plane z = the magnetic-field 
and sheet-current distributions resemble those shown 
in Figs. 2 and 3. The maximum magnetic field par- 
allel to the top surface occurs at the origin, where 
Hm(0, e,0) = (2m/j/g)i, and the first flux penetra- 
tion through the film occurs when m increases to the 
value m c = H c iy^/2. The penetrating vortex splits into 
a vortex and an antivortex, and the vortex is driven 
by the Meissner screening currents via the correspond- 
ing Lorentz force F(x,z) — <PqHm {x, e, z) to the point 
(x, y, z) = (j/o/\/2, 0, 0), while the antivortex is similarly 
driven to (x,y,z) — (— yo/V%, 0, 0). A further increase 
in the coil's current (m > m c ) will result in more nu- 
cleating vortices and antivortices and cause the devel- 
opment of vortex and antivortex flux domes centered at 
(x,y,z) — (±yo/\/2, 0, 0). Vortex and antivortex flux 
domes also could be produced by bringing a permanent 
magnet with magnetic dipole moment m = —xm close 
to the film. 

Similar effects should occur in superconducting films 
with bulk pinning. In particular, when the maximum 
magnetic field parallel to the surface, accounting for im- 
age fields, exceeds H c \ + K c , vortex and antivortex flux 
domes should be produced with properties similar to 
those described in Sees. Ill Dl and lll El In the above geom- 
etry with a small coil or a permanent magnet producing 
a magnetic moment m = -im, we expect that the main 
effect of bulk pinning will be to reduce the separation 
between the vortex and antivortex flux domes. 

In this paper we have confined our attention to the 
case in which an initially flux-free film is subjected to 
locally applied magnetic fields increasing in magnitude. 
We expect that interesting hysteretic effects, similar to 
those in a finite-width film with a geometrical barrier 
will occur when ac magnetic fields are applied. For exam- 
ple, consider a bulk-pinning-free superconducting film in 
the geometry studied in Sec.|n]but with Iq being cycled. 
As Iq is increased from zero, we expect vortex and an- 
tivortex domes to develop as predicted in Sec. Ill Cl where 
the dome boundaries a and b are determined in part by 
the condition that ^(0,6) = H c i [Eq. 1(25)1] . Suppose 
that these values are a max and b max when Iq increases 



to some maximum value I m ax, where I c \ < I max < 2J c i, 
a ma x = xl/b max > x , and b max < x — ^/yTyl- If h 
is now decreased, the condition that H x (0,e) = H cl is 
replaced by the condition that the magnetic flux under 
each flux dome remains constant; i.e., Kc[Q(a) — G(b)] 
= const. As a result, as Iq decreases, the width of each 
dome increases while its height decreases; i.e., b decreases 
(b < b max ) and a increases (a > a max ). The resulting val- 
ues of b and a can be calculated as functions of Iq using 
the relations a — x\jb and 

IoAg(a,b) — I ma xAg(a m ax,b m ax), (45) 

where [see Eqs. (|B14I) and (IB15t ] 

Ag(a,b) = Re{[g {a, y x ) - g (a, y 2 )\ 

-[g Q (b, yi )-g (b,y 2 )}}. (46) 

So long as b > 0, vortex-antivortex annihilation can- 
not occur because the sheet current flowing in the re- 
gion |x| < b still keeps vortices and antivortices apart. In 
fact, a subsequent increase of Jo would produce reversible 
changes in the width and height of the vortex and antivor- 
tex domes, provided Io < I max - However, the vortex- free 
gap of width 2b closes when Io is decreased to the value 
I ex at which b — > and a — > oo, where 

I ex = I m axAg(a max ,b max )/ Ag(co,0) (47) 

and A<?(oo,Q) = ln^/^/i)- When Iq — I ex , the sheet- 
current density becomes everywhere zero, the film ap- 
pears as if it were completely incapable of screening, and 
the magnetic field distribution (averaged over a length 
of the order of the intervortex separation) is essentially 
the same as it would be in the absence of the film. On 
the other hand, viewing the field distribution as a lin- 
ear superposition of the Meissner-state response and a 
vortex-antivortex distribution, we see that when the gap 
of width 2b closes, vortex-antivortex annihilation begins 
to occur at x = 0, and magnetic flux begins to exit from 
the vortex and antivortex domes. As Iq decreases from 
I ex to zero, the magnitude of the magnetic flux under 
each dome decreases to zero. When Iq = 0, the film 
is again flux- free, and as Iq further decreases to — I max , 
the behavior is very similar to that for increasing Iq, ex- 
cept that the roles of vortices and antivortices are inter- 
changed. 

We expect that similar but somewhat more compli- 
cated hysteretic effects will occur in the presence of bulk 
pinning. In the case of strong bulk pinning (K c 3> H cl ), 
the role of H c \ can be neglected, and the hysteretic prop- 
erties can be calculated analytically as in Ref. [25t which 
treats the response of a superconducting film to currents 
in linear wires in arrangements similar to that discussed 
in Sec. [H] Such calculations illuminate the fundamen- 
tal physics underlying the ac technique introduced by 
Claassen et alw& to determine the critical current den- 
sity j c in superconducting films. This technique employs 
a small coil, placed just above the film, carrying a sinu- 
soidal current. When the current amplitude exceeds the 
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value at which the maximum induced sheet-current den- 
sity reaches K c — j c d, a third-harmonic voltage appears 
in the coilm^i^^ A similar technique was introduced 
by Hochmuth and LorenziSi 

The effects discussed in this paper and applied to type- 
II superconducting films arc quite general and also should 
be observed in type-I superconductors. Magnetic flux 
domes consisting of intermediate-state regions contain- 
ing multiply quantized flux tubes have been observed in 
type-I strips in which the geometric barrier plays a dom- 
inant role^ It is therefore likely that separated domes 
of positive and negative magnetic flux produced in re- 
sponse to nearby current-carrying wires, coils, or per- 
manent magnets also will be observed in weak-pinning 
type-I superconducting films, foils, or plates when the 
net parallel field at the surface exceeds the bulk ther- 
modynamic critical field H c (or, when bulk pinning is 
present, H c + K c ). In type-I superconductors, however, 
we expect that the magnetic flux will enter in the form of 
the intermediate state. The analog of a vortex dome will 
be an intermediate-state region consisting either of an 
array of multiply quantized flux tubes or a meandering 
normal-superconducting domain structure carrying mag- 
netic flux up through the film, while the analog of an an- 
tivortex dome will be a similar intermediate-state region 
carrying magnetic flux down through the film. Such mag- 
netic structures should be observable by magneto-optics 
or related means by placing the magnetic-field source on 
one side of the sample and the magnetic-field detector on 
the opposite side. 
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APPENDIX A: SCREENING EFFECTS 

The primary purpose of this paper is to describe in 
detail how magnetic flux penetrates through a supercon- 
ducting film in the form of vortices and antivortices, using 
the assumption that the magnetic field below the film is 
initially negligibly small. In the following we present re- 
sults confirming that this is an excellent approximation 
when d > X and d <C y\ ■ 

Consider the experimental configuration shown in 
Fig. 1 but allow for the finite thickness d of the super- 
conducting film in the region \z\ < d/2. When the film, 
characterized by the London penetration depth A, is in 
the Meissner state, the vector potential and the mag- 
netic field throughout all space, as well as the supercur- 
rent density in the film, can be calculated as described 
in Ref. 0. The results for the x components of the mag- 
netic fields H x (0, =pd/2) at the bottom and top surfaces 
of the film, expressed in units of Ho£c(0,0) [Eq. @], are 



H x (0,-d/2) 
H Ox (0,0) 

H x (0,+d/2) 
ffox(0,0) 



Qq 



2ym 

{V2 - Vi) Jo L 2qQ cosh(Qd) + (q 2 + Q 2 ) sinh(Qd) 
Q[qcosh(Qd) + Q sinh(Qd)] 



( e -wi - e ~ qV2 )e qd/2 dq, 



2ym 

(j/2 - yi) Jo L2gQ cosh(Qd) + (q 2 + Q 2 ) sinh(Qd) 



(e 



-92/1 



-qy2\„qd/2 



dq. 



(Al) 
(A2) 



where Q = \J q 2 + A -2 . These quantities are plotted as 
the solid curves in Figs. I22landl23las functions of X/d for 
the case that d = j/i/IOOO and y 2 = 2y\. 



When A — > 0, the film screens perfectly, such that 
H x (0,-d/2) = and H x (Q,+d/2) = 2H x0 (0,d/2). 
In the opposite limit, when A — > oo, H x (0,—d/2) = 
H x0 (0,-d/2) and H x (0,d/2) = H x0 (0,d/2), where 
H x0 (0,y) = Io(y2-yx)/2ir(yi-y)(y2-y) is the magnetic 
field produced by the wires in the plane x — for y < yi 
in the film's absence. 



Equations (|A1I) and (|A2I) reduce to simpler expressions 
when d « i/i and either A < d or, if A > d, the Pearl 
length 32 A = X 2 /d obeys A <C yi- Then to good approxi- 
mation q can be set equal to zero inside the brackets and 



Q can be replaced by 1/A. The resulting integrals yield 



H x {0,-d/2) 
#0.(0,0) 

H x (0,+d/2) 
Ho* (0,0) 



2(yi + y 2 )A 
yi y 2 sinh (d/X) 
2 2{ yi +y 2 )X 
2/12/2 tanh (d/X) 



(A3) 
(A4) 



As shown by the long-dashed curves in Figs. 1221 and 12*51 

these expressions are excellent approximations, indistin- 
guishable from the solid curves, when X/d < 1 or A <C yi 
when X/d > 1. However, the long-dashed curves deviate 
significantly from the solid curves for X/d > 10, which 
is to be expected, since for the parameters used for the 
figures, A w 2/1 when X/d w 30. 

To evaluate Eqs. IjAlfl and l|A2|l for all values of the 
Pearl length^ 2 A = A 2 /c? when d <C y% and A > d, it 
is a good approximation to ignore q 2 relative to 1/A 2 
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FIG. 22: Magnetic field just below the film 
H x (0, -d/2)/H x0 {0,0) vs A/d, shown here for d = j/i/1000 
and j/2 = 2j/i as calculated from Eq. 1A11 (solid curve), 
Eq. 1A3I (long-dashed curve), and Eq. 1A5I (short-dashed 
curve) . 




FIG. 23: Magnetic field just above the film 
fl a ,(0,+d/2)/.H' !e o(0,0) vs A/d, shown here for d = j/i/1000 
and j/2 = 22/i as calculated from Eq. 1A2I (solid curve) and 
Eq. 1A4II (dashed curve). A plot of the same quantity using 
Eq. (IA6II would be indistinguishable from the solid curve. 



inside the brackets and to replace Q by 1/A but to retain 
the terms proportional to q in the denominators. This 
approximation is equivalent to the assertion that when 
d <C A, the only length that determines the screening 
properties of the film is A. This procedure yields the 
following approximate results: 

H x (0,-d/2)/H 0x (0,0) = 1-1, (AS) 
H x (0,+d/2)/H 0x (0,0) = 1 + 1, (A6) 



where 
I 



2/12/2 



(2/2 - yi) Jq 1 + 2qA 



( e -wi - e - qy2 )dq, (A7) 



^)A [eW2AEi( - 2/2A) 

-e^/ 2A Ei(- yi /2A)], 



(A8) 



and Ei(x) is the exponential integral. On the scale 
of Fig. 1221 the approximation for H x (0, —d/2) given in 



Eq. (|A5|) . shown as the short-dashed curve, is indistin- 
guishable from the solid curve for values of X/d > 3. 
Surprisingly, on the scale of Fig. 1231 the approximation 
for H x (0, d/2) given in Eq. I|A6() is indistinguishable from 
the solid curve for all values of X/d. 



APPENDIX B: COMPLEX FIELD AND 
COMPLEX POTENTIAL 

In the following we derive the complex field H(C,) and 
complex potential Q(Q satisfying the boundary condi- 
tions that the perpendicular (y) component of the mag- 
netic field in the plane of the film obeys H y (x, 0) = for 
|ar| < 6 or |x| > a and that the sheet-current density in 
the z direction obeys K z (x) = —K c for b < \x\ < a. We 
begin by considering the function 



F(C) = [W(C)-W||]MC), 

where H(() is defined in Eq. QJ, 



Wii(0 



2/i 



2/2 



and 



m = V(« 2 -c 2 )(c 2 -& 2 )- 



(Bl) 



(B2) 



(B3) 



For |£| — > oo, 0(C) - ► T*C 2 i where the upper (lower) sign 
holds in the upper (lower) half plane. Just above (below) 
the real axis, 



(x ± ie) 



±i<f>(x), \x\ < b or \x\ > a, 



<t>{x) 



b < \x\ < a, 



(B4) 



where 



\/(a 2 -x 2 )(b 2 -x 2 ), 



\x\ < b, 

4>(x) = ^ sgn(x)^/ (a 2 — x 2 )(x 2 — b 2 ), b < \x\ < a, 
-^/(x 2 -a 2 )(x 2 -b 2 ), \x\>a. 

(B5) 

As can be seen from Eqs. ||TJ and 1 132(1 . F(() is an analytic 
function of £ = x + iy except for poles at C — ±*2/i and 
£ = ±ij/2, and a branch cut along the real axis. 

Next consider the following integral around the closed 
contour C consisting of a line just above the real axis at 
£' = it + ie from it = — oo to u = +oo, the infinite circle 
at C' = Re ie from 9 = to 9 = 2tt with R -> oo, and a 
line just below the real axis at £' = it — ie from it = +oo 
to it = — oo, 



+ " du ^ + *)-^-*), (B6) 

-oo U- C 



where the integral around the infinite circle vanishes be- 
cause |F(C)| oc |i?(C)/C 2 | -> as |C| -> oo. Using the 
residue theorem, accounting for the poles of the integrand 



draft (February 6, 2008) 



14 



on the left-hand side of Eq. (|B6|I at (' = £, iyi,iy 2 , —iyi, 
and —ii/2, we obtain 



2m J_ 00 u-Q 



(B7) 



where 



Fo(Q = £(- 

27T V Si 



2/1 



2/2 



(C 2 +2/i) S2 (C 2 +2/ 2 2 ) 



;), (B8) 

si = y/W+ylW + yl), (B9) 

^2 = v /(a2+j/l)(&2 + y|). (BIO) 

From Eqs. Q, Q), and (jB2j), we find that 

H(a:±ie) -W||(a;±ie) =if y (x,Q) =Fi^(a;)/2, (Bll) 
such that Eqs. l|BTl) and l|B4|) yield 



-2iH y (x,0)/4>(x), \x\<b, 
F(u+ie)~F(u~ie) = { -iK z (x)/cf>(x), b < \x\ < a, 

-2iH y (x,0)/(j>(x), \x\ > a. 

(B12) 

However, H y (x,0) = for \x\ < b or | re j > a, and 
K z (x ) = —K c for 6 < |x| < a. Using Eqs. (|B5|) and 
(|B12() and evaluating the integrals, we obtain the follow- 
ing expressions for the complex field W(C) and complex 
potential G(0 = ff e H((')dC: 



no = -_{^r^ c^i— I 



2ttI e+yf 
±iK c /2, 



(B13) 



SiO = ^ [3o(C, 2/0 ~ So(C, 2/2)] ± ^ C, (B14) 



where 



5o(C,2/) = / rfC , 2 s 



C ,v[i ±tf£) /*] 
C' 2 + y 2 



i arctan(£/y) 

t — [a¥#(0, fc) + 2/ 2 (& 2 + y 2 )F(9, k) 

asy 

-{a 2 + y 2 )(b 2 + y 2 )II(0, ~b 2 /y 2 , k)], (B15) 

(B16) 
(BIT) 
(B18) 



8 = arcsin(C/6), and 
k = b/a, 



where E, K, and II are incomplete elliptic integrals. 

When b = 0, the following replacements can be made 
in the above expressions, 



Sl 
S-2 



2/i V a +2/i' 



2/2 V « 



2/1, 



(B19) 
(B20) 

(B21) 



and £7, X, and II can be evaluated to obtain 



g (C,y) = l arctan(C/y) + v /gg-- 7 =^ 



-arctanru/ ^ + arctanh^=2 . (B22) 



It follows from Eq. {Q that the integral / Tt(()d( 
around the great circle at |£| — > oo yields i J K z (x)dx 
along the real axis. The requirement that the film car- 
ries no net current is thus equivalent to the requirement 
that J ri(C)d(, — 0. Using Eq. I|B13(I and the property 
that for |C| — > oo, (f>(() — ► T-*C 2 , where the upper (lower) 
sign holds in the upper (lower) half plane, we obtain the 
requirement that when K c = 0, 2/i/si = 2/2/ s 2- Solving 
the latter equation, we obtain the following condition re- 
lating a and b, 



ab = x 2 = 2/12/2- 



(B23) 



Similarly, when K c > and b > 0, the requirement that 
the film carries no net current leads to the condition that 



2/1 



2/2 



V(a 2 +2/ 2 )(6 2 +2/ 2 ) V(« 2 +2/2 2 )(& 2 +2/ 2 2 ) ' ^ 

(B24) 

When K c > but b = 0, the same requirement leads to 
the condition that 



la 



n v V« 2 +2/i V^ + 



2/1 



(B25) 
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